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Abstract 

We argue that one cannot correctly calculate the elastic scattering S- 
matrix for high-energy dipole-dipole scattering, in the region where 5* is small, 
without taking fluctuations into account. The relevant fluctuations are rare 
and unimportant for general properties of inelastic collisions. We find that 
the Kovchegov equation, while giving the form of the S-matrix correctly, 
gives the exponential factor twice as large as the result which emerges when 
fluctuations are taken into account. 
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1 Introduction 



The topic of this paper is understanding the role that fluctuations play when 
unitarity corrections become important in high-energy scattering. We focus 
our attention on dipole-dipole scattering because the issues are a bit sharper 
there. We can perhaps motivate our discussion by the following observations. 
At very high energy the typical (mean) configuration of a dipole's light-cone 
wavefunction is a Color Glass Condensate [1-4], a state characterized by a 
saturation jSj momentum, Qs, and having high occupancy for all gluonic 
levels of momentum less than or equal to Qs, which depends on the rapidity, 
y, of the parent dipole. Suppose we scatter two such evolved dipoles at zero 
impact parameter, in the center of mass frame and at relative rapidity Y. 
Then, if we use just these typical configurations, of the condensate type, to 
compute the ^-matrix for the scattering, the result will be proportional to 
exp{— const Ql(Y/2)rl/a'j.}, where tq is the size of the parent dipoles serving 
as the seeds for Color Glass Condensates which collide. The Q^r"^ part of 
this formula is purely geometric; it is the number of, roughly, independent 
parts of one of the condensates when viewed on a scale Ax± ~ 1/Qs- A more 
complete discussion of condensate-condensate scattering is given in section 
3.3 with the S-matrix given in eq. (25). On the other hand the Kovchegov 
equation [7j gives an S'-matrix (see eq. (14)) which is much larger than (25), 
being proportional to exp{— const ln^((5^rQ)}. Why do eqs. (14) and (25) 
disagree and which of these answers is correct ? 

The problem is in the use of typical, or mean, configurations to estimate 
the elastic scattering S'-matrix. If one collides in the center of mass frame 
two Color Glass Condensates, then the exponent we arrived at, Q^rQ/a^, 
corresponds to the average number of gluons produced in the collision. The 
condensate description of the wavefunction is a good starting point for a 
typical (inelastic) collision, but it need not be correct for evaluating a very 
small observable like the elastic S'-matrix whose evaluation can be (and is) 
dominated by much rarer configurations. As we shall see the parts of the 
wavefunction which dominate the S'-matrix are very dependent on the frame 
in which we view the scattering. In the center of mass frame, and in the 
high-energy regime where S is very small, the S'-matrix is dominated by 
wavefunction configurations which have relatively few gluons |B,^ and which 
are not in a condensate state. In an asymmetric frame, the parent dipole 
having the lower momentum evolves into a state with few gluons while the 
higher momentum dipole evolves into a condensate, but into a condensate 
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having a lower than normal saturation momentum. 

What about the result (14) coming from the Kovchegov equation? While 
the functional form of the solution O El is fine, we find that the Kovchegov 
equation misses the correct evaluation of the exponent of by a factor of two. 
Again the problem is that the Kovchegov equation does not treat fiuctuations 
properly. The crucial ingredient in deriving the Kovchegov equation is the 
replacement of Sy\x± — z±, z± — the scattering 5— matrix for a two- 
dipole state (dipoles of sizes Xj_ — zj_ and z± — y±) by the product Sy{x± — 
Zj_)SY{z±—yj_). Such a replacement is true only in the absence of fiuctuations 
in the light-cone wavefunction of the target. Contrary to a rather widely held 
belief this replacement is not justified by large Nc behavior since, as we shall 
see below, the relevant fiuctuations concern a suppression of gluons in certain 
momentum regions of the wavefunction, an issue which is independent of Nc- 
Despite these problems, the Kovchegov equation is probably the best one can 
do with respect to a simple equation which naturally imposes unitarity on 
BFKL evolution. 

Much of what we say here has been anticipated in some early works jSlITU] 
on the study of unitarity in dipole-dipole scattering. The issue of fiuctuations 
was reasonably well understood although at that time it was not possible to 
achieve the level of analytical precision that we are able to present here. The 
detailed numerical studies done by Salam [TU] correspond to a (numerical) 
evaluation of the Balitsky equation ^T] or, equivalently, of the functional 
evolution equation for the Color Glass Condensate ^1 |H , rather than 

an evaluation of the Kovchegov equation. 

Finally, it must be admitted that we are not able to prove that the wave- 
function configurations which we suggest are dominant are in fact so. Cer- 
tainly, the configurations that we shall find are better than those chosen by 
the Kovchegov equation since our S'-matrix is much bigger. And we do give 
arguments in sections 3.6 and 4 to the end that natural modifications of our 
choice of configurations lead to a smaller S'-matrix. Nevertheless, it is diffi- 
cult to classify all the possisble configurations so as to be sure that we have 
made the right choice. In this context it should be noted that the coefficient 
of the (a^y)^ term in the exponent of our result (31) is about a factor of 1.7 
higher than the numerical evaluation in Ref. 0. We do not understand the 
resolution of the discrepancy. 

In Sec. 2 we give a rapid derivation of the Kovchegov equation, and of 
its prediction for the high-energy behavior of the S'-matrix for dipole-dipole 
scattering in the regime where S is very small. The usual Levin-Tuchin result 
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[SI El is recovered. 

In Sec. 3, we investigate dipole-dipole scattering at very high energies in 
the center of mass (CM) frame and at fixed impact parameter. In Sec. 3.1 
we review the BFKL-dipoIe picture. In Sec. 3.2, we show that using BFKL 
evolution for the dipole wavef unctions, and taking into account only average 
properties like their respective dipole number densities, leads to an ^-matrix 
which is too small, even though we work in a rapidity region where BFKL 
evolution is correct for average properties of the wavefunction. In Sec. 3.3 we 
collide two Color Glass Condenses, and again find an S'-matrix which is unre- 
alistically small. In Sec. 3.4 we describe which (rare) confingurations lead to 
the result of the Kovchegov equation in a natural asymmetric frame. In Sec. 
3.5 we describe a set of rare fluctuations which give a much larger S'-matrix 
and which we suggest dominate very high-energy scattering in the CM frame. 
In Sec. 3.6 we give (incomplete) arguments as to why the configurations we 
have chosen are optimal. 

In Sec. 4, we transform our center of mass picture to an arbitrary frame 
and describe the picture which emerges there. 

2 The Kovchegov equation 

The Kovchegov equation is probably the best "simple" equation for deal- 
ing with the onset of unitarity in high energy, but weak coupling, scattering 
in QCD. The argument for the Kovchegov equation is easily stated. Consider 
the high-energy scattering of a dipole, consisting of a quark at x\_ and an 
antiquark at on a target which may be another dipole or a more complex 
hadron or nucleus. It is convenient to view the scattering in a frame where 
the dipole is going along the negative z-axis (left-moving) and the target is 
going along the positive ^-axis (right-moving), and where almost all of the 
rapidity, Y, of the scattering is taken up by the right-moving system. (This 
is sometimes referred to as the "dipole frame".) In this frame the process 
takes the form of an elementary (unevolved) dipole of size — y± scattering 
in the field of a (highly evolved) target. If SY{x±,y±) is the ^-matrix for 
the scattering we define the amplitude Af^y by 

Afxy=l-SYix^,y±). (1) 

Now suppose we increase y by a small amount. If this increase is given to 
the target, then in order to calculate the change in Sy, it is necessary to 
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calculate how the wavefunction of the target changes with Y. This change 
can be calculated using a functional equation derived by Jalilian-Marian, 
lancu, McLerran, Weigert, Leonidov and Kovner (JIMWLK) [T| El IT^ HI] . 
However, in order to derive the Kovchegov equation it is simpler to keep 
the rapidity of the target fixed and put the small change of rapidity into 
the left-moving elementary dipole. The dipole now has a small probability 
of emitting a gluon due to this change of rapidity. If the gluon is in the 
wavefunction of the dipole at the time it scatters on the target, then what 
scatters is a quark- antiquark-gluon system which, in the large Nc limit, can 
be viewed as a system of two dipoles. (One of these dipoles is the original 
quark and the antiquark part of the gluon while the other dipole is the quark 
part of the gluon and the original antiquark.) If the gluon is not in the 
wavefunction at the time of the scattering, it can be viewed as the "virtual" 
term which decreases the probability that the quark- antiquark pair remain 
a simple dipole, compensating the probability for the two-dipole state. One 
can put this in formulae as 



■^SY{x±,y±) 



a 



271 J ^ (xx 



[Xi_ - y±} 



z±)^[z^ - y±, 



X 



s!y\x_L,z±,y_L) - SY{x_L,y±) 



(2) 



where = agNc/ir, z± is the transverse coordinate of the emitted gluon, 

(2) 

and Sy stands for the scattering of the two dipole left-moving system on 
the target. 

If the target is homogeneous on a scale large compared to x± — y±, then 
it is convenient to view Sy as a function of x± — y± and 6j_ = ^^^^^ , and Sy^ 
as a function of x± — y±, — y± and b^. Suppressing the 6^-dependence, 
eq. (2) becomes 



d_ 
dy 



Sy{xi_ -y±) 



2^ 



X 



Sy\x±- z±,z±-y±) - SY{x±-y_L) .(3) 



{x± - z±y{z± - yi_Y 



This is not yet the Kovchegov equation, rather it is part of a set of equations 
derived by Balitsky (llj , and which also follows from the functional evolution 
equation in Refs. [H |2l CH]- Eq- (3) is very difficult to use because a 
solution for the F-dependence of Sy requires knowing Sy' . One obtains the 
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Kovchegov equation if one further assumes 



Sy\x± - z_L, z± - yx) = Sy{x± - z±)Sy{z± - y±), (4) 

which is a sort of mean field approximation for the gluonic fields in the target. 
In his original papers [Zj, Kovchegov has obtained eq. (4) for the case that 
the target is a large nucleus; in that case, this equation should be a good 
approximation so long as Y is not so large as to cause S to be dominated by 
nuclear density fluctuations. However for general targets one cannot expect 
(4) to be exact, as we shall see in the next section. Nevertheless, (4) may 
be a reasonable approximation as unitarity corrections are just becoming 
important, and in any case it leads to an interesting equation. 
Specifically, using the approximation (4) in eq. (3) gives 

T^SY{x±-y±) = — d z±- (5) 

dy 2n J {x± - z±y{z± - y^y 



X 



[Sy{xi_ - z^)Sy{z± - y^) - Sy{xi_ - y±)] , 



or, equivalently. 



d ,r as f r2 ix±-y± 

-■Nxy = — d Z_L 



\2 



dy' 2nJ ^{x^-z^y{z^-y^y 

X + Ky - Ky - Kz Ky] , (6) 

which are two common forms of the Kovchegov equation [7J. Eq. (6) is the 
more useful equation when scattering is weak. In that case the quadratic 
term in (6) may be dropped and the dipole version of the BFKL equation 
results. Eq. (5) is easier to use when 5* is small, and unitarity corrections 
have become very important. 

But in the general case, when neither 5* nor A/" is small, eqs. (5) or (6) 
are difficult to deal with analytically. There are a number of good numerical 
evaluations [14-19] of the Kovchegov equation which do cover the region 
where S is neither close to zero nor close to one. Here, however, we are 
interested only in the high-energy regime where unitarity corrections are 
very important, so S is small indeed, and the solution to eq. (5) can be 
evaluated analytically. If S'y(xx — y±) is very small, then the quadratic 
term in S on the right-hand side of (5) is much smaller than the linear 
term unless either {x± — z±y or {z± — is as small as Qj'^{Y). We have 
introduced the saturation momentum, Qs{Y), which is an intrinsic scale of 
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the target (characteristic of the gluon density there), and which marks the 
scale at which a dipole scattering off the target makes the transition from 
weak (rj_ <^ to strong {r± 3> 1/Qs) interactions {r±_ is the size of 

the dipole). More precisely, it is common to define Qs{Y) by the equation 
'S'y(rj_) = 1/e when r± = 2/Qs. 

The regime of interest here corresponds to rg = \x± — y±\ ^ ^/Qg. 
Then, the right-hand side of (5) is dominated by the logarithmic regions of 
integration where one has either 

« (xx - z^'f « rl (7) 

or 

« {z^ - y^f « rl (8) 
One easily finds 

-ln^y(ro) ^-as ^— = -aM\Qlijy^. (9) 

Now, we know that (see, e.g., |1]) 

ln[Ql{Y)rl]=cas{Y-Yo) + --- (10) 
where Yq is such that QsiY^) ~ I/tq, and 01201 

MM ^ 4.883. (11) 
1 — Aq 

In the equation above, Aq is the solution to 

x(Ao) = -x'(Ao)(l-Ao), (12) 

where 

X(A)=^(1)-1^(A)-1^(1-A) (13) 



is the usual BFKL \21\ l22j eigenvalue function. The omitted terms on the 
right-hand side of (10) have only a logarithmic dependence on Y. Using (10) 
in (9) one finds il El 

5y(ro)=e-f'^^(^-^«)'^y„(ro), (14) 
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where S'y(,(ro) ~ 1. 

Because we have not kept \n(Y — Yq) terms on the right-hand side of (10), 
there are some missing hnear terms in Y in the exponent of eq. (14). Our 
focus in this paper will be primarily on the quadratic term, and its coefficient, 
in the exponent of (14). 

Eq. (14) gives the standard result in the literature. We have gone through 
such a detailed "derivation" of (14) because the main point of the present 
paper is to argue that (14) is in fact not quite right. The problem, as we 
shall see, is not with the derivation of (14) from (5), but with the mean 
field approximation contained in (4). Deep in the saturation region, where 
S is very small, fiuctuations away from average configurations are actually 
dominant for the calculation of S. 

3 Dipole-dipole scattering in the CM system 

In this section we consider dipole-dipole scattering in the center of mass 
frame. Of course, the scattering cross section cannot depend on the frame 
in which it is calculatd, however, the form that the calculation takes can be 
very different in different frames in a non-covariant gauge calculation. In the 
first three parts of this section we proceed in a qualitative manner, in order 
that the physics issues not be obscured by technical details. In the final three 
parts of this section we proceed more technically and arrive at an evaluation 
of the S'-matrix which disagrees with (14), being a factor of two smaller in 
the exponent. 

3.1 The BFKL approximation 

Suppose we scatter two dipoles, each of size ro, at a relative rapidity Y. Then 
in a frame where one of the dipoles has rapidity y and the other has rapidity 
Y — y the BFKL approximation gives |21] 



where n{ro,ri,y) is the number density of radiated dipoles of size ri in the 
wavefunction of a parent dipole of size tq, within a rapidity interval equal 
to y. Furthermore, is the lowest order dipole-dipole scattering cross 




(15) 
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section, given by 



a(o)(ri,r2) - 27ra,V| ( 1 + ln^ 1 , (16) 



where r< is the smaller of ri, r2 and r> is the larger of ri, r2. In the saddle 
point approximation to the BFKL solution, 

leading to 

g(ap-l)Y 

<7(y) 47ra^ro' , (18) 



In these equations, ap — 1 — 4Q:sln2 is the BFKL 'intercept'. What we 
are interested in here is not so much the precise formulae contained in eqs. 
(15)-(18), but rather the qualitative picture. To that end we rewrite the 
cross-section, eq. (18), as 

a{Y) = T:rlaln\r^, Tq, Y/2)C(asY) (19) 

where 



C{asY)^8^^-asN,C{S)Y (20) 

is an uninteresting and slowly varying prefactor. Eq. (19) gives the cross 
section in terms of an elementary dipole-dipole cross section, a^TiVQ, times 
the number of evolved dipoles in the wavefunctions of each of the colliding 
parent dipoles. Eq. (19) is a rehable formula so long as cr is well below ttTq. 
When a approaches hTq, unitarity corrections become important and it is to 
this topic that we next turn. 

3.2 The onset of unitarity corrections 

Let Yq be the rapidity where cr{YQ) is equal to vrrg. Then, roughly, 

^o-^-rln^, (21) 
ap — 1 a^ 
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where the neglected terms in (21) are of size In Inl/a^. Eq. (17) gives 
the average number of evolved dipoles in the parent dipole wavefunction so 
long as y < Yq, at which point saturation effects in the dipole wavefunction 
become important. (The fact that the wavefunction of a dipole does not 
receive significant saturation corrections in the region Y <^ Yq [25 follows 
from eq. (15), which if we take y small requires that the result (18) come 
completely from n{ro,r2,Y). This also shows that, when the scattering is 
seen in an asymmetric frame, unitarity corrections and saturation effects in 
the wavefunction of the evolved dipole start to be important at the same 
rapidity^, namely at Y = Yq.) 

Now consider dipole-dipole scattering in the CM frame with the two 
parent dipoles scattering at zero impact parameter and at rapidity Y < 2Yq. 
What is the S-matrix for such a collision? Of course when Y > Yq we expect 
very strong unitarity corrections to the scattering. However, it should be 
reasonable, so long as Y < 2Yo, to assume that the individual wavefunctions 
of the colliding dipoles are still given by BFKL evolution (since y = Y — y = 
Y/2 < Yq). Now Sy has the interpretation of being the probability that 
no interaction take place in the collision. It is tempting to estimate Sy by 
requiring that no evolved dipole in one wavefunction interact with any dipole 
in the other wavefunction. This gives 

In SUro) - - Co a^iro, tq, Y/2) (22) 

where the constant factor cq is not under control. When Y ~ 21^0, using (17) 
and (21) in (22) gives 

^2yo(ro)^e-^/"' (23) 

which is much smaller than what one would obtain using the result (14) from 
the Kovchegov equation along with the estimate (21). 

Eq. (23) is clearly not right. The reason for such a failure is the fact 
that the above calculation of the S-matrix in the presence of unitarity cor- 
rections, cf. eqs. (22)-(23), has included only typical configurations in the 
wavefunctions of the parent dipoles. Of course, the use of typical configura- 
tions is correct as long as the energy is not too high [Y < Yq), so that one 
can restrict oneself to the single scattering (or "single pomeron exchange") 

^This is in agreement with the discussion in Sec. 2, where we have seen that, in the 
dipole frame, the critical rapidity for strong scattering Yq is related to the saturation scale: 
Qs{Yo) ^ 1/ro, cf. eq. (10). 
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approximation, eq. (18). But at higher energies {Y > Yq), where multi- 
ple collisions are important, the wavefunction correlations among the dipoles 
which scatter simultaneously play a crucial role, and lead to a result for S 
which is very different from eqs. (22)-(23). So long as y < 2Yo, and for 
center of mass scattering, the relevant correlations are correctly described by 
BFKL evolution, but this has to be applied to the detailed dynamics of all 
the dipole configurations, and not only to average properties, like the mean 
dipole number density (17). 

The correct calculation which sums up "multiple pomeron exchanges" 
and replaces eqs. (22)-(23) in the regime where Iq < ^ < 2Yo is actually 
known. This has been originally developed within the color dipole picture 
1211 123; and recently rederived within the color glass formalism ^B] . 
This is a rather complex calculation which cannot be analytically completed, 
but has been implemented numerically by Salam jtlOj. As mentioned in the 
Introduction, the numerical results thus obtained |H] are consistent with the 
functional form in eq. (14), but with a smaller coefficient in the exponent. 
As we shall see later in this section, such a smaller coefficient emerges indeed 
when some rare configurations in the wavefunction are taken into account. 

3.3 The collision of two Color Glass Condensates 

In order to emphasize that typical wavefunction configurations can give very 
wrong answers when 5* is small, we now estimate, from typical configura- 
tions, the ^-matrix when Y ^ Yq- Then, in the CM frame, what collide are 
two Color Glass Condensates, each characterized by a saturation momentum 
Qs{Y/2). The number of gluons having k± ~ Qg in each wavefunction is 



where, as before, vq is the size of the parent dipole in each of the colliding 
condensates, ci is a constant which is not important for our purposes. Again 
interpreting Sy as the probability that no collision takes place one easily gets 



where the factor Q^Tq in (24) and (25) counts the number of independent 
regions availabe for gluons, with kj_ ~ 1/Qs, to occupy in an area vrrQ. 



Nc^c,-^Ql{Y/2)rl 



(24) 



s 




(25) 
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Alternatively, the result in eq. (25) could be obtained by using the formalism 
developed in Ref. j2ni for the scattering between two color glasses, and 
focusing on the typical configurations, of the condensate type, in both dipoles. 

When Y 2^0, Qs(>"/2) ^ (5f(>o) ^ and eq. (25) agrees with 

the previous estimate in eq. (23). Eq. (25) of course must also be wrong 
because it gives an ^-matrix which is much too small. The problem with 
(25) is the same as with (23): rare, rather than typical, configurations of 
the wavefunction dominate the evaluation of the S'-matrix in center of mass 
scattering when Y >Yq. 

3.4 Rare vs typical fluctuations in the Kovchegov equa- 
tion 

In order to better appreciate the fact that the role played in the scattering 
by the various parts of the wavefunction depends on the choice of frame we 
note that, even if we restrict ourselves to typical configurations, as in Sees. 
3.2 and 3.3, the final result for the S'-matrix can be very different if the 
calculation is performed in a different frame. For instance, let us choose the 
asymmetric frame previously used in Sec. 2, in which one of the participants 
in the collision is an elementary dipole while the other participant carries 
most of the total rapidity, and thus is highly evolved. If in this frame we 
compute Sy as the scattering between the elementary dipole and the typical 
configuration in the energetic projectile, which is a condensate, one obtains [3] 
the same result as from the Kovchegov equation, namely Eq. (14). (Indeed, 
one can view Eqs. (7) and (8) as imposing the condensate condition on 
the wavefunction of the high energy projectile in our discussion in Sec. 2.) 
Although this result is not quite right, as we shall see in the next sections, 
it is nevertheless much larger (and thus closer to the correct result) than the 
result obtained by working in the center of mass frame, Eq. (25). 

This discussion helps explain why the Kovchegov equation provides a 
quahtatively correct result although its derivation focuses on typical config- 
urations, as we have seen in Sec. 2. But this also shows that if one tries to 
interpret the result of the Kovchegov equation in a different frame, like the 
center of mass frame, this interpretation will generally correspond to some 
rare configurations, although not necessarily the optimal ones. 

Besides the "dipole frame" discussed above and in Sec. 2, where the 
relevant configurations are those of the condensate, there is another frame 
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in which the configurations pertinent to the Kovchegov equation are easily 
identified. For two parent dipoles of the same size tq, this is the frame in 
which the rapidity is equal to Y — Yq for the right-mover and to —Yq for 
the left-mover. We shall assume that the left-mover is subjected to normal 
evolution, as described by the JIMWLK equation, in the whole rapidity 
interval from —Yq to 0. Thus, at the time of scattering, this is a Color Glass 
Condensate with saturation scale Qs{Yq) = I/tq. If Y = Yq, i.e., if the right- 
mover were an elementary dipole, we would have an S-matrix of order one: 
^Yoi^o) ~ 1- The configuration which, for Y > Yq, gives the result of the 
Kovchegov equation is the one in which the evolution of the right-mover is 
suppressed in such a way that, in the whole rapidity range < y < Y — Yq, 
the corresponding wavefunction consists only of the parent dipole of size rg. 
This is a rare configuration, but has the advantage to involve only one dipole, 
and thus give a rather large contribution to S, of order one. Even though 
suppressed by the small probability of the particular configuration, that we 
shall shortly compute, this contribution remains substantially larger than 
that of a typical configuration, with N(Y — Yq) ^ 1 gluons (cf. Eq. (24)), 
which is extremely small (of the same order as shown in Eq. (25)). 

Let A{rQ, Y — Yq) denote the probability of the rare configuration of in- 
terest. This is the same as the survival probability of the parent dipole after 
a BFKL evolution over a rapidity interval Y — Yq, and can be computed with 
the methods of Refs. [HI HDl I2S1 1211 123 I2i • This probability decreases with 
increasing Y because of gluon radiation, and the corresponding rate is the 
same as the virtual term in Eq. (6): 

^ A(ro, y) = -aJ^ ^ A{rQ, y) (26) 

oy J 2tc {x± - zx_)\z^ - y±)^ 

with x± and yj_ labelling the quark and antiquark parts of the dipole vq. This 
equation should be integrated from y = up to y = Y — Yq. Note however 
that the integral is ill-defined because of poles in the integrand at z± = x± 
and z± = yj_. These divergences refiect the fact that one cannot avoid the 
emission of dipoles of arbitrarily small size. 

The situation is similar to that encountered in the discussion of the 
Kovchegov equation in the regime where the S-matrix is small (cf. Sec. 
2). Like in that physical cutoff, which depends on y, exists also for the 

problem at hand. Let us introduce this cutoff as a minimal dipole size p{y); 
that is, the integral in (26) should be restricted to {x± — z±)^ > p^iy) and 
(zj_ — ?/_l)^ > p^iy). Physically this means that, at rapidity y, we suppress 
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Figure 1: The configuration retained by the Kovchegov equation in the frame 
in which the left mover has rapidity —Yq. 

the radiation of dipoles with sizes larger than p{y), but smaller dipoles are 
allowed^. Each of these small radiated dipoles becomes the seed of a normal 
BFKL evolution going from the intermediate rapidity y (at which the dipole 
was emitted) down to y = 0. The evolution is such that new dipoles are 
produced, and the typical size of such "children" dipoles is increasing with 
decreasing y. This is so because the line ln{r^{y)/rl) = —easy, where c is 
given in (11), is a line of constant amplitude for BFKL evolution P71 128j . 

We see that, instead of being truly a single dipole configuration, the 
configuration that we are constructing also contains small dipoles, of size 
< p{y) at the "moment" y of their emission, which then grow up with further 
decreasing y, and can become as large as p^(0) = p^(?/)e'^"*^ at the collision 
"time" y = 0. Still, this complicated configuration behaves, during scattering, 
as the idealized single dipole configuration, provided all the additional dipoles 
which are present in the wavefunction at ?/ = 0, and hence also in the interval 

■'Note indeed that the dipoles which are included in the integral in cq. (26) are those 
whose radiation is forbidden. 
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< y < Y — Yq, are much smaller than the parent dipole Tq. (Indeed, such 
small dipoles undergo essentually no scattering, so the overall S'-matrix for 
this configuration remains very close to 5*1^, (ro) ~ 1.) The condition that 
p(0) ^ ro implies an upper limit on p{y), namely p'^{y) ^ rge"'^"''^, which 
can also be written as 

ln{p\y)/rl) ~ -ccx.y, (27) 

since the integral on the right-hand side of (26) is dominated by the loga- 
rithmic regions of integration where p{y) ^ r ^ tq. Here, r is the transverse 
size of any of the emitted dipoles, that is, either \x± — z±\, or \y± — z±\. 
Thus, the dominant contribution to eq. (26) can be evaluated as: 



— lnA{ro,y) ~ -a^ / — = asln{p^{y)/rl) = -ca^y, 

oy Jp2(y) r 



(28) 



which after integration over y yields: 

A{ro, Y-Yo) ~ e-i"'(^-^°)'. (29) 

Of course, the evolution towards larger dipoles, with sizes r ^ tq, must be 
forbidden as well. But such large dipoles give only a small contribution to 
the integral in eq. (26), because of the rapid fall off of the integrand at large 
transverse separations l-z^ — ^ Tq (with b± — {x± + y±)/2). Specifically, if 
we denote this contribution as A' (this is an additional factor which multiplies 
the previous contribution in eq. (29)), then A' is estimated as: 

dy ^ 2 ° 2 

After integration over y, this gives a new contribution, of order ds{Y — Yq), 
to the exponent of eq. (29). This contribution is subleading in the regime 
of interest here, and, in any case, it goes beyond the general accuracy of the 
present calculation. Because of that, dipoles with sizes r > ro will be ignored 
in what follows. 

The 5*— matrix element associated with the particular configuration of 
interest is finally computed ets Sy — A{ro, Y — Yq)Syo coincides, as 
anticipated, with the result of the Kovchegov equation, Eq. (14). 

In the next subsection and in Sec. 4 we shall see that a larger S'-matrix 
can be obtained by choosing different rare configurations, which are also 
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frame dependent. To that end, it will be useful to have a graphical rep- 
resentation of the regions of evolution in the (ln(/c^rQ), F) plane, with fc^ 
being the momentum conjugate to r±. For the particular configuration that 
we have discussed above, the corresponding plot is displayed in Fig. 1. The 
shaded triangle in this plot represents the kinematical domain into which 
dipole emission is forbidden, and which gives the dominant contribution to 
A, eq. (29). Note that, up to a factor dg, the exponent in eq. (29) is the 
same as the area of this shaded triangle. 

3.5 The dominant configuration in the center of mass 
frame 

We now turn to the task of finding those configurations that dominate very 
high energy dipole-dipole scattering in the CM frame. Since the typical 
configurations that we have previously considered, in Sees. 3.2 and 3.3, have 
given an S'-matrix which is too small (cf. eqs. (23) and (25)), it is natural 
to search for configurations which are more rare in the wavefunction but 
which lead to a larger S'-matrix. The problem with the previous calculations 
is that the typical configurations contain too many gluons at the time of 
collision, thus giving a contribution to the S'-matrix which is extremely small. 
This suggests that the dominant configurations for evaluating S at very high 
energies should contain less than the mean number of gluons 8J. 

In view of the discussion in Sec. 3.4, the general strategy for finding the 
optimal configurations should be rather clear: Loosely speaking, the final 
goal should be to minimize the number of gluons by suppresing the evolution, 
but at the same time maximize the probability of the resulting configuration. 
The latter is largest (i.e., of order one) for the configurations subjected to 
normal evolution, but is rapidly decreasing for the configurations in which 
the evolution has been suppressed. Thus, the best results should be obtained 
by enforcing the minimal suppression which still give rise to a S'-matrix of 
order one. 

Unfortunately, since we are not able to analyze all the possible configura- 
tions, we cannot transform these considerations into a precise mathematical 
criterion for the search of the optimal configurations. What we can do, how- 
ever, is to use the experience with the previous example in Sec. 3.4 in order 
to guess what should be the dominant configurations. With these config- 
urations, we shall then evaluate the S'-matrix, and find a result which is 
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Figure 2: The optimal configuration in the center-of-mass frame. 

indeed larger than in eq. (14). Finally, in later sections, we shall give argu- 
ments as to why wc believe we have found the most important wavefunction 
configurations for evaluating S in the high energy limit. 

Consider a zero impact parameter collision in the center-of-mass frame 
at rapidity Y > Yq. (As usual, Yq denotes the critical value for the onset 
of unitarity corrections, cf. eq. (21).) For this problem, we shall require 
that the wavefunction of the right-moving dipole consist only of the parent 
dipole, of size ro, in the rapidity interval lo/2 < y < Y/2, with a similar 
requirement on the wavefunction of the left-moving dipole in the interval 
-Y/2 < y < -Yo/2. In the rapidity intervals -Yq/2 < y < for the left- 
moving system, and, respectively, < y < lo/2 for the right- moving one, we 
allow normal, BFKL, evolution of the wavefunctions. 

Of course, say in the Yq/2 < y <Y/2 interval for the right-moving system, 
we cannot require that all evolution in rapidity be absent. What is required 
is that evolution which does exist create only very small dipoles, close to 
either the quark or antiquark of the parent dipole, so that in the interval 
Yq/2 < y < Y/2 the system have no more than one dipole of size Xtq or 
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larger, with A a constant of order 1. In order to guarantee that this be the 
case, it is necessary to suppress the creation of dipoles much smaller than 
To at rapidities y > Yo/2; otherwise, the dipoles emitted at intermediate 
rapidities could evolve into dipoles of size tq, or larger, at rapidity Yo/2. 
The constraints are similar to those we have just discussed in the previous 
subsection and are illustrated in Fig. 2. Gluon emission from the two parent 
dipoles, as part of the evolution which forms the left and right-moving states 
which scatter on each other, is forbidden if the gluon has k± and y values 
lying in the shaded triangles*^ of Fig. 2. (For the moment ignore the small 
unshaded square in the upper triangle.) The line 

ln{klrl) = ca{y-Yo/2), (30) 

and a similar line for the lower triangle, is determined, as we have just seen, 
by the requirement that gluons lying to the right of that line, for any y > 
Yo/2, cannot evolve through normal BFKL evolution to give gluons having 
A;j_ < 1/ro at rapidity Yo/2. The constant c is again the same as in (10) and 
so the hue (30) is a hne of constant amphtude in terms of BFKL evolution. 
Thus evolution which crosses this line has a small probability. 

Still as in Sec. 3.4, we shall denote by A{ro, (Y — Yo)/2) the probability 
that the parent dipole, of size tq and rapidity Y/2, not give rise to any emis- 
sions in the upper triangle of Fig. 2. The ^'-matrix is then given by a factor 
A{ro, {Y — Yo)/2) for each of the parent dipoles partaking in the collision, 
times the (partial) S'-matrix for the scattering of two dipoles separated by a 
rapidity gap Yo/2 — {—Yo/2) — Yo, and which are subjected both to normal 
evolution. The latter is, of course, S'yo(ro) ~ 1. After also using eq. (29) for 
A{ro, {Y — Yo)/2), one finally obtains: 

5y(ro)~e-t^^(^-^o)^5^„(ro) (31) 

a result which is much larger than (25), and also significantly larger than the 
result (14) coming from the Kovchegov equation. 

ThTis, one can avoid a very small S— matrix, as given in (25), by not for- 
getting the rare parts of the wavefunction which can dominate the 5*— matrix 
when S is small. These rare configurations consist of unusually small num- 
bers of gluons being present in the wavefunction. 

^As discussed in Sec. 3.4, the evolution towards larger dipoles, with sizes r > ro 
(or In(fc^rQ) < 0), is also forbidden, but this can be neglected in the calculation of the 
probability A{ro, {Y - Fo)/2). 
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It should be emphasized that (31) holds only when the parent dipoles 
have a fixed size, Tq in our case. For instance, if one were to scatter heavy 
onia states, eq. (31) would not emerge. In the case of heavy onia there is a 
wavefunction giving the probability of having a parent dipole of a particular 
size, r. Because of that, the yS-matrix for the scattering of two onia at very 
high energies would be dominated by parent dipole sizes r of order l/Qs{Y). 

3.6 Sampling other configurations 

In the last section we found a particular configuration that leads to an 
5"— matrix which is much larger than that given by the Kovchcgov equa- 
tion. While this is sufiicient to prove that fiuctuations must be important, 
and thus that the result in eq. (4) cannot be exact, it is, perhaps, not clear 
that we have found the optimal configurations for dipole-dipole scattering in 
the CM frame. In this section we give arguments to the end that (31) is the 
correct answer and that the configurations we have focused on are indeed the 
dominant ones. 

The first change one may consider in the previous calculation would be 
to modify the values of the intermediate rapidities ±Yo/2 which separate 
between suppressed and normal evolution. So, let us replace Yq ^ Yi, with 
Yi Yq, in Fig. 2 and the related calculations. There are two possibilities 
— Yi < Yq, or Yi > Yq — , but as we argue now they both lead to a S~ 
matrix which is much smaller than that in eq. (31). i) If Yi < Yq, the 
probability A{ro, {Y — Yi)/2) for the suppressed configuration is even smaller 
than before, while the partial S'-matrix Sy^ for the collision between two 
normally evolved systems is still of order one (since Yi < Yq corresponds 
to weak scattering). Thus, this situation is clearly less favourable, ii) If 
Yi > Yq, the probability A^Tq, (Y — Yi)/2) is larger, but this gain is more 
than compensated by the strong decrease in the partial ^'-matrix Sy^^- To 
see this, note that for the calculation of Sy^ we are either in the situation 
described in Sec. 3.2 (if Yi < 2Yo), or in that of Sec. 3.3 (if Yi > 2yo), 
and, as we have seen, both these situations give only tiny contributions to 
S. If we focus on | lnS'y|, for definiteness, then | InSyJ grows exponentially 
with the difference Yi — Yq (cf. eqs. (22) or (25)), while the corresponding 
decrease in \\nA\ is only quadratic. Clearly, the overall (S'-matrix decreases 
rapidly with increasing Yi — Yq. 

The other major change that we could make in our calculation in the last 
section would be to allow some limited emission of gluons into the forbidden 
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Figure 3: The shaded area is the domain of suppressed evolution for a parent 
dipole within the small unshaded rectangle in Fig. 2. 

triangles of Fig. 2. Recall that the area of such a triangle determines the 
penalty factor A associated with suppressing the evolution (cf. the remark 
at the end of Sec. 3.4). Then let us calculate the change in our previous 
result which occurs if we allow emission into the small square in the upper 
triangle of Fig. 2. Take the area of the square to be l/cts, so that now one 
gets a factor A which is a factor of e larger than the previous one: 

A' = exAiro,iY-Yo)/2), (32) 

with A{ro, ^2^° ) given by eq. (29). We have less suppression than before 
because we have not forbidden the parent dipole to emit into the small square. 
However, now we must forbid the newly created (small) dipole from evolving 
further and again creating a strong suppression like that given in (25). In 
Fig. 3, the shaded triangle represents the region into which one must forbid 
the additional dipole from emitting gluons^. If y and k± are the coordinates 
of the gluon emitted by the parent dipole (i.e., the central coordinates of the 
small unshaded square in Fig. 2), and if yi is defined by (see Fig. 3): 

\n{klrl) = casiyi-Yo/2), (33) 

^Once again, we consider only that forbidden area which gives the dominant contribu- 
tion to the probabihty A; see eq. (35). 
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then the shaded area in Fig. 3 — the area into which emission from the 
secondary dipole must now be forbidden — is 

a = ^ca^iy -yi)'^ . (34) 
The probabihty of not emitting into this area is: 

In order for this suppression not to be larger than the gain from being allowed 
to have the secondary emission, cf. eq. (32), one must require that a^a < 1. 
This implies the following condition on the location of the small square in 
Fig. 2 : 

^ < 2c, (36) 

where, for more clarity, we have denoted the rightmost borderline of the 
suppressed area in Fig. 3 as (cf. eq. (30)): In ((5^(z/)^o) = '^^iv ~ ^o/2). 

Eq. (36) shows that, for a given y, the point ln(A;^rQ) lies within one unit 
of the boundary. In {jQ'^ {y)r'^ . But since the boundary is ambiguous at this 
level of precision — e.g., the right hand side of eq. (30) is specified only up 
to corrective terms of order one — , we see that it does not seem possible to 
relax the assumption we made in deriving eq. (31), namely that there must 
be no emission into the triangle regions of Fig. 2. 



In {Q\y)rl) - In {klrl) 



4 The picture in a general frame 

In this section we show how the result (31) comes about in an arbitrary 
frame. In order to describe the frame choice and the scattering picture, it is 
useful to refer to Fig. 4. We scatter a left-moving parent dipole of size ri 
and rapidity —Y2 on a right-moving parent dipole having size ro and rapidity 
Y — Y2. Recall that Iq is the value of the rapidity difference between two such 
dipolcs where the S-matrix begins to be significantly different from one. We 
also suppose that F2 < ~ ^0), for later convenience. 

Clearly one cannot allow both dipoles ro and ri to have any significant 
amount of normal BFKL evolution, or else we will find a strong suppression 
of the type given by (25). Since the left-moving dipole has the smaller rapid- 
ity it is natural (and easier) to suppress its evolution, and we will describe 
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^ln(k2r2)=cas(y-Yi-Yo) 



^In (k2r2)=cas(Yi-y) 



In (k2 r 2) 



Figure 4: Choosing the optimal configuration in some generic frame. 



that suppression in a moment. As for the dipole ro, we suppress evolution 
over its Y — Y2 — {Yi +Yo) highest units of rapidity exactly as we did in Sec. 

3, with the region of suppression given by the upper shaded triangle of Fig. 

4. In the lowest Yq + Yi units of rapidity the parent dipole undergoes normal 
evolution. We will determine Yi later, by maximizing the ^'-matrix. The 
unshaded triangle, whose rapidity values go from to denotes the satu- 
ration region for the right-mover, that is, the region where the parent dipole 
ro has evolved into a Color Glass Condensate (so the evolution is non-linear 
within that particular region). Let Qs = QsiYo + ^1) be the saturation scale 
which characterizes this condensate. Since obtained after a normal evolution 
through 10 + ^1 units of rapidity, this is of the form: 

Ql{Yo + Y,) = g^e'=^(^»+^^), (37) 

where Qg intrinsic scale of the right-moving system, proportional to 

l/rg and, possibly, powers of Og. By the very definition of Yq, we have 
Q1{Yq) = 1/rf, so the saturation momentum (37) can be rewritten as: 

Ql ^ QHYq + Y,) = ie^^^S (38) 
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where the dependence upon the unknown scale Qq has disappeared. 

Turning now to the parent dipole ri, we require that no additional dipoles 
be created, through gluon emission, which would have a strong interaction 
with the right-moving Color Glass Condensate. Clearly, the dipoles which 
would have such strong interactions are those which, at the time of scattering 
(i.e., a,t y — 0), would be of size r ^ 1/Qs- This means that, at all the 
intermediate rapidities within the range —Y2 < y < 0, we have to suppress 
the emission of those dipoles which, after a normal evolution over the last \y\ 
units of rapidities, could become of size l/Qs, or larger. By referring to the 
similar discussion in Sec. 3.4, we see that the maximal size p{y) allowed for 
a dipole emitted at rapidity y is constrained by 

P^y) ^ T^e-'^"^!^' = r^e-'^"^^^^-^), (39) 

Qs 

where we have also used eq. (38) together with the fact that \y\ = —y. In 
terms of the conjugate momentum variable k±, this gives the lower shaded 
region in Fig. 4 as the region into which radiation is forbidden. 
We are now in a position to estimate the S'-matrix as 

5y(ro,ri) = Aii{ro,Y-Yo-Yi-Y2)SY,+Y,{ro,ri)AL{n,Y2), (40) 

where and Al are the suppression factors from the no emission require- 
ments for the two dipoles, and 5yo+y^(ro, ri) is the 5*— matrix for scattering 
an elementary dipole of size ri on a dipole tq which has evolved into a Color 
Glass Condensate characterized by the saturation momentum (38). 

To compute S, one can rely on the Kovchegov equation, or, more exactly, 
on its approximate form in eq. (9), valid deeply at saturation. This is correct 
for the present purpose since, as discussed at the beginning of Sec. 3.4, the 
Kovchegov equation describes correctly the scattering between an elementary 
dipole (here, the dipole ri) and the typical configurations in a Color Glass 
Condensate (here, the dipole tq which has been allowed to carry out its 
normal evolution over the rapidity interval Yq + Yi). By using eq. (14), one 
obtains: 

SYo+Yriro,n) ^ e-§^?^i'5y„(ro). (41) 

We note in passing that the exponent in (38) is, except for a factor of ds, the 
area of the unshaded triangle in Fig. 4, and where, as usual, the exponents 
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of Aji and Al are given in terms of the area of the upper and lower shaded 
regions of Fig. 4. Thus 

An = e-t^'(^-^^-^^-^°)', (42) 

and, similarly, 

By combining eqs. (40)-(43), one obtains: 

- exp {-^ al [{Y -Y,-Y,- Y,f + {Y, + Y^f] } Sy,,{ro) ■ (44) 

The S'-matrix given by (44) refers to a particular set of configurations 
of the wavcfunctions, characterized by one parameter: the rapidity Yi which 
describes the amount of evolution in the right-moving system. (Recall that 
the other rapidity parameter in eq. (44), namely Y2, specifies the frame of 
reference.) The actual S'-matrix should be determined by the value of Yi 
which maximizes the right hand side of eq. (44), or, equivalently, which 
minimizes the exponent there. This condition yields: 

Yi^liY- Yo) - Y,. (45) 

This is non- negative, because of our initial assumption Y2 < ^(Y — Yq). 

Finally, Sy is obtained by evaluating eq. (44) with the optimal value for 
Yi, eq. (45) : 

SY{ro,n) = e-t^»(^-^«)'Sy„(ro). (46) 

(The right hand side depends upon ri via the "critical" rapidity Yq.) This 
result is independent of Y2 (i.e., upon the choice of a frame), and is exactly 
the same as the corresponding result in the CM frame, eq. (31). 

In fact, the picture of the evolution for the optimal configuration in the 
CM frame, as shown previously in Fig. 2, can be obtained from the general 
picture in Fig. 4 by first choosing Y2 in such a way that Yi = — this requires 
Y2 = (Y — Yq)/2 (cf. eq. (45)) — , and then performing a supplimentary boost 
by a rapidity amount AY = —Yq/2 (so that Y2 — > 5^/2, as in Fig. 2). The 
latter operation changes only the position of on the rapidity axis, but not 
also the picture of the evolution^. 

^Indeed, if Yi = 0, there is no non-linear evolution involved, but only linear, BFKL, 
evolution, which is invariant under a small boost. 
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In particular, it is instructive to visualize the evolution in the "dipole 
frame", i.e., the rest frame of the left-mover {Y2 = 0), in which both the 
optimal configuration (cf. Fig. 4) and the configuration retained by the 
Kovchegov equation are explicitly known. These configurations are illus- 
trated in Figs. 5. a and b, respectively. As usual, a shaded area represents 
a kinematical domain into which dipole emission is forbidden, while an un- 
shaded triangle represents the saturation region. Also, up to a factor of fig, 
the total area of the (shaded and unshaded) triangles represents the exponent 
of the ^-matrix associated with the respective configuration. As manifest 
on these figures, the area of the empty triangle in Fig. 5.b is twice as much 
as the sum of the areas of the two small triangles in Fig. 5. a. 

Note also that, for both configurations in Fig. 5, the right-mover ends up 
as a Color Glass Condensate at y = 0. But while the configuration retained 
by the Kovchegov equation (cf. Fig. 5.b) corresponds to normal evolution 
over the whole rapidity range — BFKL evolution from y — Y down to 
y — Y — Yq, and then non-linear, JIMWLK, evolution over the remaining 
Y — Yq units of rapidity — , in the optimal configuration (cf Fig. 5. a), the 
normal evolution is allowed only over the lower {Y + Yq)/2 units of rapidity. 
As a result, the saturation scale which characterizes the condensate in the 
optimal configuration, namely, QsHY + ^o)/2), is much smaller than the 
corresponding scale QsiY) for the normal evolution. 
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(b) 

Figure 5: The optimal configuration (above) and the configuration retained 
by the Kovchegov equation (below) in the frame in which 12 = 0. 
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